Tangled magnetic fields and CMBR signal from reionization 

epoch 

Rajesh Gopafl and Shiv K. Sethfl 

Raman Research Institute, Bangalore 560080, India 

Abstract 

We compute the secondary CMBR anisotropy signal from the reionization of the universe in 
the presence of tangled magnetic fields. We consider the tangled-magnetic-field-induced scalar, 
vector, and tensor modes for our analysis. The most interesting signal for £ 100 arise from 
tensor perturbations. In particular, we show that the enhancement observed by WMAP in the 
TE cross-correlation signal for £ 10 could be explained by tensor TE cross-correlation from 
tangled magnetic fields generated during the inflationary epoch for magnetic field strength Bq ~ 
4.5 x 10 -9 G and magnetic field power spectrum spectral index n ~ —2.9. Alternatively, a mixture 
of tensor mode signal with primordial scalar modes give weaker bounds on the value of the optical 
depth to the reionization surface, T re ion : Treion = 0.11 ± 0.02. This analysis can also be translated 
to a limit on magnetic field strength of ~ 5 x 10~ 9 G for wave-numbers ^ 0.05 Mpc -1 . 
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I. INTRODUCTION 



Coherent magnetic field of micro-gauss strength are observed in galaxies and clusters 



of galaxies 



:Q 



for a recent review see e.g. [3|). Observational evidence exist for even 



larger scale magnetic fields Q. The origin of these observed magnetic fields however is not 
well understood. The observed magnetic fields could have arisen from dynamo amplification 




rich originated from various 



lid, m 



. Alternatively 



of small seed ( ^ 1CT 20 G) magnetic fields (see e. 
astrophysical processes in the early universe 
the magnetic fields of nano-Gauss strength could have originated from some early universe 
process like electroweak phase transition or during inflation (e.g. Q, Q, see Q, Q 
for reviews). In this scenario, the observed micro-gauss magnetic fields then result from 
adiabatic compression of this primordial magnetic field. 

The existence of primordial magnetic fields of nano-Gauss stren gth can influence the large 



scale structure formation in the universe 



;0, 0, 0, 0,0, 



Also these magnetic 

he.ds eo U .d ,eave observe signatures in the CMBR anisotropies ( M , Q, Q, Q, Q, 

3, 0)- 

In recent years, the study of CMBR anisotropies has proved to be the best probe of 
the theo.es of strnctore formation in the We.se (see , g . a for a recent review). The 
simplest model of scalar, adiabatic perturbations, generated during the inflationary era, 
appear to be in good agreement with both the CMBR anisotropy measurements and the 
distribution of matter at the present epoch (see e.g. 0], 0|). Tensor perturbations could 
have been sourced by primordial gravitational waves during the inflationary epoch. There is 
no definitive evidence of the existence of tensor perturbations in the CMBR anisotropy data; 
the WMAP experiment, from tem pera ture anisotropy data, obtained upper limits on the 
amplitude of tensor perturbations [30(. Vector perturbations are generally not considered 
in the standard analysis as the primordial vector perturbations would have decayed by 
the epoch of recombination in the absence of a source. An indisputable signal of vector 
and tensor modes is that unlike scalar modes these perturbations generate B-type CMBR 
polarization anisotropies (see e.g. |37| and references therein). At present, only upper limits 
exist on this polarization mode [32j. However, the on-going CMBR probe WMAP and the 
upcoming experiment Planck surveyor have the capability of unravelling the effects of vector 
and tensor perturbations. 
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Recent WMAP results suggest that the universe underwent an epoch of re-ionization at 
z ~ 15; in particular WMAP analysis concluded that the optical depth to the last reion- 



ization surface is r reion = 0.17 ± 0.04 



33|; which means that nearly 20% of CMBR photons 



re-scattered during the period of reionization. The secondary anisotropies generated during 
this re-scattering leave interesting signatures especially in CMBR polarization anisotropies 
(see e.g. as i s evidenced by the recent WMAP results 2^1 . 

Primordial magnetic fields source all three kinds of perturbations. In this paper we study 
the secondary CMBR anisotropies, generated during the epoch of reionization, from vector, 
tensor, and scalar modes, in the presence of primordial tangled magnetic fields. Recently, 
Lewis (2?! computed fully-numerically CMBR vector and tensor temperature and polar- 
ization anisotropies in the pre sence of magnetic fields including the effects of reionization. 
Seshadri and Subramanian 29} calculated the secondary temperature anisotropies from vec- 
tor modes owing to reionization. Our approach is to compute the secondary temperature and 
polarization anisotropies semi-analytically by identifying the main sources of anisotropies in 
each case; we compute the anisotropies by using the formalism of [37]. We also compute 
the tensor primary signal to compare with the already existing analytical results for tensor 
anisotropies [3| . 

In the next section, we set up the preliminaries by discussing the models for primordial 
magnetic fields and the process of reionization. In §3, §4, and §5, we consider vector, tensor, 
and scalar modes. In §6 the detectability of the signal is discussed. In §7, we present and 
summarize our conclusions. While presenting numerical results in this paper, we use the 

m H, Q; 



currently-favoured FRW model: spatially flat with VL m = 0.3 and fi A = 0.7 (|3J|, |3j, |40j) 
with Q b h 2 = 0.024 (Q, £|) and h = 0.7 ( 42|). 



II. PRIMORDIAL MAGNETIC FIELDS, REIONIZATION, AND CMBR 
ANISOTROPIES 

Assuming that the tangled magnetic fields are generated by some process in the early 
universe, e.g. during inflationary epoch, magnetic fields at large scales ( ^.O.lMpc) are 



not affected appreciably 
recombination universe ( 



)V different processes in either the pre-recombination or the post- 



H, [l 



In this regime, the magnetic field decays as 1/a 
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from the expansion of the universe. This allows us to express: B(x, 77) = B(x)/a ; here x 



3 



is the comoving coordinate. We further assume tangled magnetic fields, B, present in the 
early universe, to be an isotropic, homogeneous, and Gaussian random process. This allows 



one to write, in Fourier space (see .e.g. 



my- 



(A(q)S*(k)) = S 3 D (q - k) fa - k^/k 2 ) M(k) (1) 

Here M{k) is the magnetic field power spectrum and k — |k| is the comoving wavenumber. 
We assume a power-law magnetic field power spectrum here: M(k) = Ak n . We consider the 
range of scales between /c min taken to be zero here and small scale cut-off at k = /c max ; /c max is 
determined by the effects of damping by radiative viscosity before recombination. Following 
Jedamzik et al. 0, ^ ^ 60 Mpc" 1 (B Q /(3 x 1(T 9 G); B is the RMS of magnetic field 
fluctuations at the present epoch. A can be calculated by fixing the value of the RMS of 
the magnetic field, B , smoothed at a given scale, k c . Using a sharp A;-space filter, we get, 

7r 2 (3 + n) a 

A ~ k (3 + n) B o (2) 

We take k c = IMpc" 1 throughout this paper. For n ~ —3, the spectral indices of interest 
in this paper, the RMS filtered at any scale has weak dependence on the scale of filtering. 

Recent WMAP observations showed that the universe might have got ionized at redshifts 
z c± 15. However the details of the ionization history of the universe during the reionization 
era are still unknown; for instance the universe might have got reionized at z — 15 and 
remained fully ionized till the present or the universe might have got partially reionized 
with ionized fraction x e £ 0.3 at z ~ 30 and became fully ionized for z ^ 10. Both these 
ionization histories are compatible with the WMAP results [3^. Given this lack of knowledge 
we model the reionization history by assuming the following visibility function, which gives 
the normalized probability that the photon last scattered between epoch rj and 77 + drj, to 
model the period of reionization: 

9(V> Vo) = rexp(-r) = — ^ Treion )) exp \_(^ _ r/rcion ) 2 /Ar/ r 2 cion j (3) 

V7rA?7 rc i on J 

Here r(r/, r] Q ) = J^ Q n e a t dt is the optical depth from Thompson scattering; r reion is the optical 
depth to the epoch of reionization; for compatibility with WMAP results, we use r re i on = 0.17 
throughout. Rreion and A?7 re i on are the epoch of reionization and the width of reionization 
phase, respectively; we take ry reion corresponding to -2 reion = 15 and A77 reion = 0.25r7 reion . 
Notice that the visibility function is normalized to r reion for r rcion << 1. 



III. CMBR ANISOTROPIES FROM VECTOR MODES 



From a given wave number k of vector perturbations, the contribution to CMBR tem- 
perature and polarization anisotropics to a given angular mode £ can be expressed as (see 
e.g. 0): 

= ~^Io° dV eM-r)rP v miHVo - V)} (5) 
B ^ (//n - /, ' ) - -VQ T d V exp(-r)iP v ( V )(3J[k( Vo - V )] (6) 



(2£+l) Jo 

Here and V are the line-of-sight components of the vortical component of the baryon 
velocity and the vector metric perturbation. P v (r]) = l/10[6y 2 _ V^G)^] and the Bessel 
unctions, ji, €£ and that give radial projection for a given mode are given in reference 
37| . The evolution of vector metric perturbations, Vi(k, 77) is determined from Einstein's 
equations (e.g. 37 1, 2^|): 

a k 

-k 2 Vi = IQirGa 2 "£( Pj + Pi )K - ^) (8) 
j 

Here Si, the source of vector perturbations, is determined by primordial tangled magnetic 
field in this paper. The index j correspond to baryonic, photons and dark matter vortical 
component of velocities. For tangled magnetic fields, the vortical velocity component of 
the dark matter doesn't couple to the source of vector perturbations to linear order and 
Qi = vj — Vi decays as 1/a for dark matter (see e.g. |26(); and hence the dark matter 
contribution can be dropped from the Einstein's equations. The photons couple to baryons 
through Thompson scattering. In the pre-recombination epoch, the photons are tightly 
coupled to the baryons as the time scale of Thompson scattering is short as compared to 
the expansion rate; besides the photon density is comparable to baryon density at the epoch 
of recombination. In the reionized models we consider here, neither the photons are tightly 
coupled to baryons nor are they dynamically important. Therefore photon contribution can 
also be dropped from Eq. (jSJ). Eq. (jSJ) then simplifies to: 

- k 2 Vi = lQTrGa 2 p B n v B (9) 
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with Q V B = (i>£ — Vi). The quantity of interest is the angular power spectrum of the CMBR 
anisotropies which is obtained from squaring Eqs. (@]), (J5J), and (JSj), taking ensemble average, 
and integrating over all k: 

(10) 

This expression is valid for both vector and tensor perturbations; for scalar perturbation the 
pref actor is 2/tt. 

For primordial magnetic field, the sources Si(k, rj) of vector perturbation (Eq. (JJJ)) is the 
vortical component of the Lorentz force: 

S^r,) = -^xF.T.[B(x)x(VxB(x))] = S t Qc)\ (11) 

It can be checked that this Newtonian expression for Si is the same as the more rigorously 
defined Uf in Appendix A (Eq. (J2ED). 



7T 



dkk 2 







T,E,Bt 



(Mo) 



2£ + l 



A. Temperature anisotropies from vector modes 

As seen from Eq. (@J), there are three sources of temperature anisotropies. The 
most important contribution comes from vorticity Q V B . For the reionized models, using 
Eqs. (J7J) and (jHJ), it can be expressed as: 

QUKV) = ^ 1 (12) 
ap b0 

Here pbo is the baryon density at the present epoch. The other major contribution is from 
temperature quadrupole 6^ 2 . For reionized models, the quadrupole at the epoch of reion- 
ization is dominated by the free-streaming of the dipole from the last scattering surface (see 
discussion below, Eq. (115)1 ). This contribution is generally small but in this case can be 
comparable to the vorticity effects at small values of t. This is owing to the fact that the 
vorticity is decaying and therefore during reionization epoch its contribution is smaller as 
compared to the epoch of recombination. The quadrupole term on the other hand gets its 
contribution from the vorticity computed at the epoch of recombination (Eq. (|15)l). This, 
as we shall discuss below, is not the case for scalar and tensor anisotropies, as the dominant 
source of anisotropy is either constant (metric perturbations for tensor perturbations) or 
is increasing (compressional velocity mode for scalar perturbation) as the universe evolves. 
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The third source of temperature anisotropics is metric vector perturbation V; this term can 
be comparable to the other terms only at super-horizon scales. We drop this term in this 
paper. 

In Figure 1 we show the secondary temperature anisotropics generated during the epoch 
of reionization from vector modes. It is seen that the quadrupole term has significant 
contribution only for i ^ 20. The dominant contribution at larger t is from the vorticity 
during reionization. The vorticity source contribution can be approximated as: 

®h(Vo,k) S(k) y 2 .(ii) r/r )} (13) 

(c\ n i -i \ — TTT ~ Jt FWO ^reionJJTVeion U J J 
(21 +1) 47Tp b0 Tfreion 

for t < 20 and 

®TfiVo> k ) S{k) 1 r]jg{r] -£/k,r] ) fW 
(2£+l) ~ k 

oTTPbO ?7reion * *• 

for I £j 50. The temperature angular power spectrum from vorticity increases roughly as l 2A 
for i £ 50, with th e sig nal reaching a value roughly 0.3 /ik at i ~ 10 4 . This is in agreement 
with the results of 1291. 



B. Polarization anisotropies from vector modes 

The main source of the polarization anisotropies is the temperature quadrupole 0^ 2 - O ne 
contribution to the temperature quadrupole at the epoch of reionization is from the free- 
streaming of the dipole from the last scattering surface. The dipole at the last scattering 
surface can be obtained from the tight-coupling solutions to the temperature anisotropies 
. The quadrupole from the free-streaming of the dipole at the epoch of recombination 

is: 

©T2(M) = 5tt v B (k,r] rec )j!2 U \k(ri -r? rcc )) (15) 

Here 77 rec corresponds to the epoch of recombination. As is the case for scalar perturbation- 
induced polarization in the reionized model (e.g. (3^) this quadrupole doesn't suffer the 
suppression as the quadrupole prior to the epoch of recombination when the photons and 
baryons are tightly coupled. The structure of anisotropies generated by the quadrupole 
is determined by j% \krj) around the epoch of reionization. This typically gives a peak 
in anisotropies at I ~ 2r] /r] re i on . This source dominates the contribution to polarization 
anisotropies for £ £ 10. Another contribution to the temperature quadrupole at the epoch 
of reionization comes from the secondary temperature anisotropies generated at the epoch 
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of reionization. The approximate value of this quadrupole can be got from retaining the 
first term in Eq. (J3J: 

= £ drj' gfa r)')ttUv')j2 1] Hv ~ v')} (16) 

This contribution is generically smaller than the first contribution. Firstly, this depends on 
the vorticity evaluated close to the epoch of reionization as opposed to the first contribution 
which is proportional to the vorticity at the epoch of recombination. As the vorticity decays 
as a -1 / 2 in the matter-dominated era (Eq. (fT2"J) ). the latter contribution is suppressed by 
nearly a factor a 100 in the angular power spectrum. Second, as only a small fraction of 
photons re-scatter (nearly 20%), this contribution is further suppressed by a factor of T r 2 eion . 
However, this contribution is not suppressed at small angular scales and, therefore, might 
dominate the polarization anisotropies at large values of i. 

In Figure 2, we show the E and B polarization angular power spectrum from the sources 
given by Eqs. (fTHjl and (fTT)|) . As discussed above, the secondary polarization anisotropies are 
dominated by the quadrupole generated by free-streaming of dipole at the last scattering 

the -B-mode signal is larger than the .E-mode 



surface. As expected for vector modes 
signal; the signal strength reaches ~ 10 _3 /xk at i ~ 10 in both cases. This dominates the 
primary signal for t ^ 10 as also seen in the numerical results of Lewis [2^ . The contribution 
from the quadrupole generated at the epoch of reionization is seen to be completely sub- 
dominant 



IV. CMBR ANISOTROPIES FROM TENSOR MODES 

The energy-momentum tensor for magnetic fields has a non- vanishing traceless, transverse 
component which sources the corresponding tensor metric perturbation. This in turn affects 
the propagation of radiation from the last scattering surface to the present and hence gets 
manifested as additional anisotropies. In this section we calculate the effect of reionization 
on the resultant anisotropies. For the temperature anisotropies, we study this effect, by 
calculating the power spectra separately for the standard recombination (no-reionization) 
and reionized scenario whereas for the polarization anisotropies we compute the secondary 
anisotropies by using the visibility function given by Eq (j3J). 



S 



A. Tensor temperature anisotropics 



The line-of-sight integral solution for temperature anisotropics, for tensor perturbations 
is given by (j^): 

5|^2> = jTV-*[Tpm - huPw* - ,)] (17) 

Here, P T (f]) = l/10[O^ 2 — y/EQ^] * s the tensor polarization source and h is the gravita- 
tional wave contribution whose evolution is detailed in Appendix B. The polarization source 
is modulated by the visibility function and hence is localized to the last-scattering surface. 
In the tight coupling limit before recombination, P T ~ —h/(3f) for a more detailed 

derivation of P T in the tight- coupling regime see Appendix B. In the post-recombination 
epoch, P T is determined by the free-streaming of quadrupole generated at the last scattering 
surface. However, the visibility function is very small at epochs prior to reionization. There- 
fore the main contribution of this term comes only from epochs prior to recombination. The 
gravitational wave source on the other hand being modulated by the cumulative visibility 
exp(— r) contributes at all epochs. As a result, the P T contributes negligibly to temperature 
anisotropies at all multipoles for the case of standard recombination. In the reionized model, 
this term gets additional contribution from epochs close to reionization redshift but contin- 
ues to be sub-dominant to the other term. We have also checked this numerically. Hence 
we can neglect the first term in the above solution and using the matter-dominated solution 
for h (Appendix B) we arrive at the following expression for the angular power spectrum: 

Here, x = kr], x = krj and Xd = kr] rec . The above expression is evaluated numerically for 
the two different ionization histories: standard recombination with and without reionization 
which are essentially characterized by the different behaviour of the cumulative visibility 
exp(— r). The temperature power spectra are shown in Figure 3. As seen in the figure, the 
temperature power spectrum in both cases shows similar behaviour. The power is nearly 
flat upto I ~ 100 after which the amplitude falls rapidly. This behaviour is identical to 
that obtained for primordial gravitational waves. This is expected because, the tensor met- 
ric perturbation is sourced by the magnetic field only upto the neutrino-decoupling epoch 
thereby imprinting an initial nearly scale-invariant spectrum after which the evolution is 
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source-free. The effect of reionization is to reduce the cumulative visibility between recom- 
bination (z ~ 1100) and reionization (z ~ 15) epochs. This is why the signal is suppressed 
for the reionized model. Approximate analytic expressions to primary Cj T were derived in 
( 3)- However these give the correct qualitative behaviour Cf T oc £ ' 2 only for I <, 100. This 
is because in their analytic results, the lower limit for the time-integral is taken to be zero 
in Eq. (|18|) whereas the correct lower limit is r) rec since the cumulative visibility is zero for 
V ~ ?7rec- We have not neglected this lower limit in our numerical calculation and hence we 
obtain the damping behaviour for i £ 100; Lewis [27[ also obtains this damping behaviour. 
Our results are in reasonable agreement with the results of [^J in the entire range of l\ these 
results also agree to within factors with the results of (3] for I <, 75 when the different con- 
vention we use for defining B is taken into account. Our results are quantitatively accurate 
to better than 10% for the lower multipoles I ^ 75 but begin to differ appreciably from the 
results of numerical studies for larger i or in the damping regime 4^. This is because we 
have not treated the transition regime from radiation-dominated to matter-dominated for 
the gravitational wave evolution accurately. As described in Appendix B, we have assumed 
instantaneous transition. This however does not affect the qualitative description of modes 
whose wavelength is greater than the transition-width k £ rf^ which in turn corresponds to 
multipoles £ £ 300. 

B. Polarization anisotropies from tensor modes 

The line-of-sight solution for the E and B-mode polarization is given as: 

^Mh^. = _V6 jT d V feM-r)P T PHKVo - v)} (19) 

= -V6 jT drjf eM-r)P T eJ[k(Vo ~ v)} (20) 

Here, j3f and ef are the tensor polarization radial functions as given in Ref . 0] • The tensor 
polarization source P T (ri) in this case will contribute significantly to the above integral only 
close to the reionization epoch. There are two contributions to the polarization at r] VC i on : 
one due to the quadrupole generated at the reionization surface and the other due to the 
free-streaming primary quadrupole. However, as in the case of vector perturbations, the 
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free-streaming primary quadrupole will give the dominant contribution. We thus have, 

P T (k,v) = ^©? 2 (M) = ~\ I" drjhA 22) [k( Vo - V )]exp(-r) (21) 

To simplify the calculations we make the following approximation. Since the visibility func- 
tion is strongly peaked at ?7 re i n, we take P T outside the integral by evaluating it at the 
visibility peak r] re i OI1 . We have verified that this approximation works extremely well for the 
lower multipoles where the power is significant. We thus get the following expressions for 
the polarization angular power spectra: 
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fi f r -i 2 / fVO \ ^ 

Cj B = - / dkk 2 U 2 T (k) P T (r ?rcion ) / d-nte- T t3j[k{ m - rj)} (22) 

Cj E = - / dkk 2 U 2 T (k) \P T (Vreion)} / d^e^ ef [k( Vo - T,)} (23) 

71 J L \Jrirec J 

As seen in the above expressions, the polarization power spectrum is modulated by the 
visibility function itself instead of the cumulative visibility in the case of temperature power 
spectrum. As a result, both E as well as B mode anisotropies peak close to the multipole 
corresponding to the horizon scale at reionization. Physically this can be understood as 
follows: the modes which are super-horizon at reionization experience negligible integrated 
Sachs- Wolfe effect before ^ re ion and hence very small polarization is generated for such modes. 
Maximum polarization is generated for modes that just enter the horizon at ^reion- For 
sub-horizon modes, the amplitude of the gravitational wave falls and then sets itself into 
oscillations which is reflected as a drop in power for higher multipoles. 

The polarization power spectra are shown in Figure 4. As seen in the figure, the E- 
mode power peaks at I ~ 8 whereas the B-mode power peaks at t ~ 7. The corresponding 
signal strengths at the peaks are ~ 0.2fiK in both cases. As expected, the E-mode power is 
marginally greater than the B-mode power mainly because of the slightly different behaviour 
of the radial projection factors (37 1. The primary anisotropies for both the polarization 
modes is sub-dominant on these scales. This enhancement in the net (primary+secondary) 
signal was also seen in the numerical calculations of Ref. 27 1. 

We also show the primary CMBR polarization anisotropies from tensor modes in Figure 4. 
For computing these anisotropies we use the tight-coupling quadrupole, P T {rj) as derived 
in Appendix B (Eq. (|45jl ). The primary power spectra are also computed from Eqs. ()22|) 
and (}2*3*j) with lower limit of the time integral replaced by zero. Our results are in agreement 
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with the numerical results of |27[ when we take into account the fact that we use different 
value of 7]*/r)i n (Eq. (JJ3J)): we use r}*/rji n = 10 18 , which gives the epoch of generation of 
the tangled magnetic field close to inflationary epoch. While presenting numerical results, 
Lewis |2jj uses r]*/r]- m = 10 6 , which puts the epoch of generation of magnetic field close 
to the epoch of electro-weak phase transition. Therefore our signal is roughly an order of 
magnitude larger than the results of 0]. 

In Figure 5 we show the expected TE cross-correlation from tensor modes, computed using 
Eqs. (fT7|). (pOjl. and (JUJ), including the effect of reionization. The effect of reionization is 
seen as the peak in the TE cross-correlation for £ £ 10. The signal is dominated by the 
primary signal for large multipoles. Note that the TE cross-correlation is positive in the 
entire range £ £ 150 as was also pointed out by Ref. [26j] for the primary tensor TE cross- 
correlation (for details of sign of TE cross-correlation for various modes see 0]). In the 
next section we compare this signal with the WMAP observation of TE cross-correlation. 

V. SECONDARY CMBR ANISOTROPIES FROM SCALAR MODES 

In addition to the vortical component of the velocity field, the tangled magnetic fields 
also generate compressional velocity fields which seed density perturbations. These density 
perturbations have interesting consequences for the formation of structures in the universe 
([3, 13, 0], Ql' 01' Si)- The compressional velocity field also give rise to secondary 
anisotropies during the epoch of reionization. We compute this anisotropy here. The line- 
of-sight solution to the temperature anisotropies from these velocity perturbations is: 

= rdve-Tvl(k, V )jriHvo-v)] (24) 



2£+l jo 

Here vt is the line-of-sight component of the compressional velocity field, defined in 

Ref. [37]. The growing mode of compressional velocity can be expressed as (0], 13 ) : 

v s h (k, rj) = (F.T.[B(x)x(VxB(x))]) = v° h (k) V (25) 

Here p m o is the matter density (baryons and the cold dark matter) at the present epoch. 
The compressional velocity field, unlike the vortical mode, has a growing mode. Also unlike 
the vortical mode (Eq . (|12|l). the compressional mode of baryonic velocity couples to the 
dark matter (0, 0). 



12 



In Figure 6 we show the angular power spectrum of the secondary temperature 
anisotropies generated by the compressional velocity mode. The signal has a peak at roughly 
the angular scale that corresponds to the width of the visibility function during reionization 



(for detailed discussion see e.g. 44j). The amplitude of this secondary anisotropy is several 
orders of magnitude smaller than the observed temperature anisotropies and it is unlikely 
that this signal could be detected. 

VI. DETECTABILITY 

It follows from Figure 1 to 5, that the most important signal at small multipoles arises 
from tensor polarization anisotropies. In particular, the yet-undetected 5-mode signal holds 



the promise o 
authors (e.g. 



unravelling the presence of primordial magnetic fields, as also noted by other 
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]) In Figure 4, we show the expected errors on the detection of polarization 
signal from the future CMBR mission, Planck surveyor. The expected la error, valid for 
t < 100, is (e.g. Q, G) : 

AC -(pirW) (c ' + ^ (26) 

For Planck surveyor, / s k y — 1 and w ~ 1.7 x 10 16 for one-year integration. In Figure 4, we 
use the primordial tensor 5-mode signal for calculating the expected la error from Eq. (|26|). 
Figure 4 shows that the signal from magnetic fields with strength £ 3 x 10~ 9 G is detectable 
by this future mission. However, it is likely that, except for the B mode signal, the magnetic 
field signal will be buried in a larger signal. However, owing to the non-Gaussianity of the 
mag netic field signal it might still be possible to extract this component of the signal (e.g. 

H). 

In Figure 5 we show the TE cross-correlation signal from tensor modes along with the 
expected signal from primordial scalar modes with_7v eion = 0.17, which is in good agreement 
with the WMAP data of TE cross-correlation 33]. It could be asked if the TE cross- 
correlation observed by WMAP for I £ 100 could be explained as the tensor signal. From 
Figure 5 it is seen that the tensor signal at small multipoles is roughly a factor of 5 smaller 
than the scalar signal. And therefore, as the power spectrum from tangled magnetic fields 
ex Bq, much of the enhancement observed in the TE cross-correlation for I £j 10 could be 
explainable in terms of the tensor signal from primordial magnetic field for B ~ 4.5 x 10 -9 G. 
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We quantify this notion by computing the x 2 f° r ^ < 15 for both the best fit model from 
WMAP and the tensor model with B ~ 4.5 x 10~ 9 G against the detected WMAP signal 



33j | : the x 2 P er degree of freedom in the two cases is ~ 1.7 and ~ 1.8, respectively. 



Therefore the enhancement can entirely be interpreted in terms of the secondary signal from 
primordial magnetic fields. 

A more realistic possibility is that both primordial scalar and tensor modes gave compa- 
rable contribution to the observed signal. As the strength of both these signals for £ ^ 15 is 
roughly oc T r 2 eion (for details of secondary scalar signal see e.g. |38(), and assuming that there 
is roughly equal contribution from both, the inferred value of r re i n from the analysis of the 
signal could be smaller by a factor of a/2. To quantify this statement, we did a x 2 test to 
estimate r re ion by adding the tensor signal with Bq — 4.5 x 10~ 9 G and the primordial scalar 
signal with the best-fit cosmological parameters from WMAP. From this analysis we obtain 
r reion ~ 0.11 ± 0.02 (la) with a determined by 5x 2 = 1. A possible test of this hypothesis is 
non-gaussianity of the signal at small multipoles, as the magnetic-field-sourced tensor signal 
is not Gaussian. 

The tensor signal (primary plus secondary) could be appreciable for £ ;$ 100. In the 
range 15 £ £ £ 100, the tensor and primordial scalar signals are nearly independent of the 
value of T reion . While the primordial scalar TE signal anti-correlates for £ ^ 40, the tensor 
signal shows positive cross-correlation in the range £ £ 100, as seen in Figure 5. The present 
WMAP data shows tentative detection of TE anti-correlation for £ £ 100 j^. From x 2 
analysis in the range 15 ^ £ ^ 100, we notice that the tensor signal alone is a poor fit to 
the data (x 2 per degree of freedom of ~ 2.1 as opposed to a value of 1.6 for the primordial 
scalar model). However a sum of these two signals with B - 4.5 x 10~ 9 G is a reasonable 
fit, as it is dominated by the primordial scalar signal. 

It should be noted that for B Q ~ 4.5 x 10~ 9 G, the tensor temperature signal is comparable 
to the primordial scalar signal (Figure 3). WMAP analysis obtained an upper limit of ~ 0.7 
on the ratio of tensor to scalar signal ( l 30]). While this limit is rather weak, a more detailed 
analysis of the temperature signal including the effect of tensor mode signal sourced by 
primordial magnetic fields might give independent constraints on the strength of primordial 
magnetic fields. 

In our x 2 analysis we use only the diagonal components of the Fisher matrix. However, 
owing to incomplete sky, the signal is correlated, especially for small multipoles, across 
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neighboring multipoles. However, a more comprehensive analysis taking into this correlation 
is likely to yield similar conclusions for the reasons stated above. 

Our conclusions are not too sensitive to the value of small scale cut-off k max or the 
scale of the filter k c used to define the normalization (Eq. ({2)1 ) for magnetic field power 
spectrum index n = —2.9 we use throughout the paper. For k max = k c = 0.05 Mpc -1 , the 
foregoing discussion related to tensor mode anisotropics would be valid for Bq ~ 5 x 10 -9 G. 
Therefore, the results for TE cross-correlation from tensor perturbations can be interpreted 
to put bounds on magnetic fields for only large scales k £ 0.05 Mpc -1 . 

The strongest bound on primordial magnetic fields arises from tensor perturbations in 
the pre-recombination era 48]. These bounds are weakest for nearly scale invariant (n ~ 
—3) magnetic fields power spectrum (Eq. (33) of Q|) and largely motivated the choice of 
the power spectral index we consider here. For n = —2.9, the bound obtained by 3] is 
considerably weaker than B Q ~ 4.5 x 10~ 9 G, the values of interest to us in this paper. 
Vector modes might leave observable signature in the temperature and polarization signal 
for £ ^ 2000; the current observations give weak bound of Bq ^ 8 x 10~ 9 G j^J. Tangled 
magnetic-field-sourced primary scalar temperature signal gives even weaker bounds [50]. 
More recently, Chen et al. j49( obtained, from WMAP data analysis, a limit of ^ 10~ 8 G 
on the primordial magnetic field strength for nearly scale invariant spectra we consider 
here; |49j consider vector mode temperature signal in their analysis and study possible non- 
Gaussianity in the WMAP data. Another strong constraint on large scale tangled magnetic 
fields comes from Faraday rotation of high redshift radio sources (see e.g H); this constraint 
is also weaker than the value of magnetic field required to explain the enhancement of the 
TE cross-correlation signal as seen by WMAP 19j . Therefore, the value of B required to 
give appreciable contribution to the TE signal is well within the upper limits on Bq from 
other considerations. 

It should be noted that the entire foregoing discussion on the tangled-magnetic-field 
tensor signal can be mapped to primordial tensor modes. The reason for this assertion is that 
magnetic fields source tensor modes only prior to the epoch of neutrino decoupling, and the 
subsequent evolution is source free, which is similar to the primordial tensor modes which are 
generated only during the inflationary epoch and evolve without sources at subsequent times. 
Therefore, an analysis similar to ours could be used to put constraints on the relative strength 
of the tensor to scalar mode contribution (for a fixed scale) and the tensor spectral index 
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of the primordial modes. The main observational difference between such an interpretation 
and the one give here is that tensor signal sourced by magnetic fields will not obey Gaussian 
statistics as opposed to the primordial tensor modes. 



VII. SUMMARY AND CONCLUSIONS 

We have computed the secondary anisotropies from the reionization of the universe in 
the presence of tangled primordial magnetic fields. Throughout our analysis we use the 
nearly scale invariant magnetic field power spectrum with n = —2.9. For vector modes, 
we compute the secondary temperature and E and B mode polarization auto-correlation 
signal. For scalar modes, the results for secondary temperature angular power spectrum from 
compressional velocity modes are presented. For tensor modes, in addition to the secondary 
temperature and polarization angular power spectra, we compute the TE cross-correlation 
signal and compare it with the existing WMAP data; we also recompute the primary signal 
for tensor modes. Whenever possible we compare our results with the results existing in 
the literature. In particular, Lewis [^tJ recently computed fully-numerically the vector and 
tensor primary and secondary temperature and polarization power spectra. We compare our 
semi-analytic results with this analysis and find good agreement. Seshadri and Subramanian 
[29! ] computed the secondary temperature anisotropies from vector modes. Our results are 
in good agreement with their conclusion. Mack et al. j3| computed primary signal from 
vector and tensor modes using the formalism we adopt in this paper. Our results are in 
disagreement with their results for I ^ 75, and we have given reasons for our disagreement 
in the discussion above. In addition to comparison with existing literature, we also give new 
results for secondary TE cross-correlation from tensor modes and secondary temperature 
angular power spectrum from scalar modes. 

We discuss below the details of expected signal from each of the perturbation mode: 
Vector modes: The secondary temperature and polarization signals from the vector modes 
is shown are Figure 1 and 2. The secondary temperature signal increases oc £ 2A for I ^ 50 
and reaches a value ~ 0.1 (yuk) 2 for i ~ 10 4 , in agreement with the analysis of 2^|. For small 
t the signal is very small 10~ 4 (yuk) 2 ) and for large i the secondary signal is smaller than 
the primary signal (e.g. S^]) and therefore it is unlikely that the signature of reionization 
could be detected in the vector-mode temperature anisotropies. The polarization signal, 
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shown in Figure 2, is sourced by the free-streaming of dipole at the epoch of recombination. 
This signal dominates the primary signal for £ £ 10, but is several orders of magnitude 
smaller than the expected signal from tensor modes. 

Scalar modes: We only compute the secondary temperature anisotropies from compres- 
sional velocity modes in this case. As seen in Figure 6, this contribution is several orders 
of magnitude smaller than the already-detected primary signal and therefore its effects are 
unlikely to be detectable. 

Tensor modes: As seen from Figures 4 and 5, the most interesting CMBR anisotropy 
signal for £ £ 100 is from these modes. The secondary S-mode signal from tensor modes is 
detectable by future CMBR mission Planck surveyor for B ~ 3 x 10~ 9 G . The tensor TE 
cross correlation from primordial magnetic fields can explain the observed enhancement of 
the observed signal for £ 10 by WMAP for Bq ~ 4.5 x 10~ 9 G if the primordial magnetic 
fields are generated during the epoch of inflation. Assuming that tensor modes make a 
significant contribution to the observed enhancement, the bounds on the optical depth to 
the surface of reionization, r re i on are weaker by roughly a factor of V2. This hypothesis can 
be borne/ruled out by testing the Gaussianity of the signal for £ £ 10 
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Appendix A 



In this section, we briefly discuss the terminology and present the complete expressions 
for the vector and tensor power spectra U v (k) and U T (k) ([26]). The energy momentum 
tensor for magnetic fields for a single Fourier mode is a convolution of different Fourier 
modes and is given by: 



r«(k) 



d 3 q 



BiiqjBjQz - q) - ^£ m (q)£ m (k - q) 



(27) 



The energy-momentum tensor has non-vanishing scalar, vector, and tensor components. The 
vector and tensor components, in Fourier space, are defined as: 



n 



V 



P h T 

1 vp n q- L pq 



(28) 
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Here = dij — k; L kj. The vector and tensor anisotropic stress are then defined as the 
two-point correlations of the above components as: 

(nT(k)nr(k')> = 2|n^(fc)| 2 ( 5(k + k') (so) 

(n5(k)n5(k'))=4|n( T )(A;)| 2 5(k + k') (31) 



By evaluating the above correlations as also given in [26|, we can arrive at the following 
approximate expression for the power spectra for n < —3/2. 

^ - ^irrr 3 (32) 

Here, A is the normalization of the magnetic power spectrum given in Eq. (j2J) 



Appendix B 

Gravitational waves correspond to transverse, traceless perturbations to the metric: 8gij = 
2a 2 (r])hij with ha = kihij = 0. Since hy is a stochastic variable we can define its power 
spectrum as: 

(%(k,7 7 )/ iii (k',r 7 )) = 4|/ i (A;,7 7 )| 2 5(k + k') (34) 
The evolution of hij then follows from the tensor Einstein equation (see e.g. j^) , 

h + 2-h + k 2 h = 87cGS(k, V ) (35) 

(X 

The source on the RHS is the tensor anisotropic stress of the plasma which is defined 
as: S(k,i]) = U T (k)/a 2 (Eq. (|3"Tf ). We assume that the primordial magnetic fields are 
generated by some mechanism at a very early epoch r\i n . It was recently shown by Lewis [^J 
that after the neutrino decoupling epoch 77* the neutrino start free-streaming and develop 
significant anisotropic stress which cancel the anisotropic stress of the primordial magnetic 
fields to the leading order for super-horizon modes, resulting in negligible net anisotropic 
stress in the plasma. We can thus assume that for 77 ^> 77*, S^k, 77) = and for 77 <C 77*, 
S(k,i]) = U T (k)/a 2 where U T (k) is the magnetic tensor anisotropic stress as defined in 
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Eq. (|3ip. We now derive the solutions to Eq. (J33j) in various regimes. The evolution of the 
scale factor a{rj) is given by the Friedmann equation: 

a 2 = H*(n m a + tt^ + tt u + tt A a*) (36) 

Here, fl m ,y,u,A are the fractional densities in matter, radiation,neutrinos and cosmologi- 
cal constant respectively. Approximate solutions in the radiation-dominated and matter- 
dominated epoch are 0(77) = 2y / a ^"" "-^ an d a{rf) = respectively. Using the above 
form for the scale-factor we can rewrite Eq. (|3*5|) for r] in < r] < 77* as: 

h + -h + k 2 h = — 2 y —L (37) 

V Py V 

Here, i? 7 = f2 7 /(f2 7 + f2 v ) ~ 0.6. p 7 is the CMBR energy density. Eq. (}3*Tj) can be solved 

exactly using the Green's function technique to give 



3: 



m ^ = 3B,nr ( t) r. sin[fc(>; - t/)] 
For super- horizon modes fcf|<C 1, the above form can be simplified to give: 

h(M) „ r ^ = ffl f (39 ) 

For 77 3> r/*, the evolution of ft, is given by the homogeneous solutions in the radiation and 
matter-dominated regimes: 

h ra d(k,v) = Ajo{kv) (40) 
h mat (k,r ] )=A 2 3 -^- (41) 

The coefficients A\ and A 2 are determined by matching the super-horizon solutions at the 
two transitions 77* and T] eq . We thus get 

A 2 = SA, = 9R - aT(k) h, f » ) (42) 
Thus, the full expression for the matter-dominated solution can be written as: 



Ut) .?tflJi)*) (43) 

Pi VW V 

This solution is used for solving tensor temperature and polarization primary and secondary 
anisotropies. Few assumptions have been made in deriving the above expression. Firstly, 
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the transition between radiation dominated to matter-dominated region has been assumed 
to be instantaneous. This however does not affect the evolution of modes with wave-length 
greater than the width of transition ki] ccl £ 1. Moreover, only super-horizon solutions have 
been used to match the solutions for h at different transitions. These simplifications however 
do not affect the results quotes for small multipoles as discussed in the main section. 



Tight-coupling tensor quadrupole 

In the tight-coupling regime, z £ 1100, to lowest order in mean-free path, we have P T = 
—h/(3r) |26j]. We however use the expression accurate to the second order in mean-free path 
as is done for the scalar modes in [51]. Using the Boltzmann equation for the evolution of 
tensor modes we get the following equation for P T (k, rf) in the tight-coupling limit: 

P + —fP = -— (44) 
10 10 V 1 

The lowest order solution to this equation is obtained by neglecting the P in the equation, 

which gives, P = —h/(3f). The above equation however can be solved exactly to give: 

P(rj) = r drfhe-^^-*™ (45) 
Jo 

We use the standard recombination history for computing r. 
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FIG. 1: The secondary temperature angular power spectrum from vector modes is shown. The solid 
and the dashed lines correspond to the contribution from vorticity and the total signal, respectively 
(see text for details). The power spectrum is plotted for Bq = 3 x 1CP 9 and n = —2.9 (Eq. @). 
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FIG. 2: The secondary polarization angular power spectrum from vector modes is shown. The 
solid and the dot-dashed lines correspond, respectively, to the B and E mode contribution from 
the free-streaming quadrupole (Eq. H15j0 . The dotted and dashed curves B and E mode signals 
that arises from the source term given by Eq. (|16|) , The power spectra are plotted for Bq = 3 x 10 -9 
and n = -2.9 (Eq. Q). 
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FIG. 3: The contribution of tensor modes to the temperature power spectrum is shown. The 
solid and dashed lines give, respectively, the power spectra without and with reionization. The 
thick solid line, shown here for comparison, correspond to the temperature power spectrum from 
scalar modes, for the best-fit parameters from WMAP ( 
B = 3 x 1(T 9 , n = -2.9 (Eq.Q), and r/*/r? in = 10 18 (Eq. gfl) 



3(J|). The power spectra are plotted for 
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FIG. 4: The tensor secondary and primary polarization power spectra are shown along with the 
expected signals from primordial scalar and tensor modes. The dot-dot-dot dashed and thin dotted 
lines correspond to the secondary E and B mode power spectra, respectively. The dot-dashed and 
dashed lines give the primary E and B mode power spectra. The two top solid lines, shown here 
for comparison, correspond to the E and B mode power spectra from primordial scalar modes, for 



the best-fit parameters from WMAP ( 30J). For B mode signal we assume the ratio of tensor to 
scalar quadrupole T/S = 0.7 and the tensor spectral index rit = 0. The bottom solid lines shows 
the .B-mode signal expected from gravitational lensing. The thick dashed line shows the 1-a errors 
expected from the future CMBR experiment Planck surveyor for one year of integration (Eq. (|26[)). 
The power spectra are plotted for B = 3 x 1(T 9 , n = -2.9 (Eq.©), and rj^/vin = 10 18 (Eq. (E310 . 
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FIG. 5: The tensor TE cross-correlation power spectra are shown along with the expected signal 
from primordial scalar modes. The dot-dashed shows the TE cross-correlation (secondary plus 
primary) from tangled magnetic fields. The thick solid line, shown here for comparison, correspond 
to the (absolute value of) TE cross-correlation power spectrum from primordial scalar modes, for 
the best-fit parameters from WMAP ([s(|). The power spectrum is plotted for Bq = 3 x 1CP 9 , 
n = -2.9 (Eq.©), and r)*/r) in = 10 18 (Eq. (@3J)). 
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FIG. 6: The secondary temperature power spectrum from scalar mode perturbations, seeded by 
tangled magnetic fields, is shown. The power spectrum is plotted for Bq = 3 x 10 -9 and n = —2.9 
(Eq.©). 
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